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We review the latest information that is available about the parton distributions of the 
■ , proton, paying particular attention to the determination of the gluon. We briefly describe 

the various processes that have been advocated to be a measure of the gluon. We discuss 

(— ! 

the importance of the gluon to the description of the structure function F% at small x, 
with emphasis on the lnl/x resummations. 
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1. Parton distributions 



Perturbative QCD is remarkably successful in describing the broad sweep of hard scattering 
processes involving the proton. A vital common ingredient is a universal set of parton distri- 
butions, fi(x,Q 2 ), which allow all of these reactions to be calculated in terms of basic QCD 
subprocesses at the partonic level. fi(x, Q 2 ) is the probability of finding parton i (where i may 
be a quark, antiquark or gluon) within the proton carrying a fraction x of its momentum when 
probed by a particle with virtuality Q 2 . 

The classic way to probe the partonic structure of the proton is deep-inelastic lepton-nucleon 
scattering, where the lepton may be an electron, a muon or a neutrino. At high energy the 
differential cross-section for, say, deep-inelastic electron-proton scattering (ep — > eX) has the 
form 

with Q 2 = —q 2 , the Bjorken x variable x = Q 2 /2p.q and y = Q 2 /xs, where p and q are 
the 4-momenta of the proton and virtual exchanged photon respectively, yfs is the centre-of- 
mass energy of the electron-proton collision. It is easy to see that the momentum fraction x 
is the same as the Bjorken x. Since the struck quark acquires 4-momentum xp + q we have 
{xp + q) 2 = m 2 . Thus x = Q 2 /2p.q in the infinite momentum frame where masses may be 
disregarded. 

The relation between the observable structure function F 2 and the parton densities is, to 
0(as), of the form 




F 2 (x,Q 2 )=xy2e 2 / £ lf q {£,Q 
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where the partonic subprocesses are shown in Fig. 1. The 0(«s) QCD subprocesses shown in 
(b,c) have initial state collinear singularities, which are factored off into the parton densities 
causing them to "run" (i.e. to depend on Q 2 ) leaving well-behaved known coefficient functions 
C q and C g . Due to this renormalisation, the absolute values of the parton densities are not 
calculable in perturbative QCD. Rather QCD determines the Q 2 dependence (or so-called 
scaling violations). It is given by the DGLAP evolution equations [0] which have the form 

dfj(x,Q 2 ) _ ^ f dx' D fx 
In 

where the splitting functions 



= Pii®fi, (3) 



P l3 =a s pV + a 2 3 PV + ... (4) 

So far, the leading order (LO), P. fj , and next-to-leading order (NLO), Pjj \ terms have been 
calculated. 
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Figure 1: Partonic subprocesses which contribute to deep-inelastic scattering: (a) the lowest-order 
diagram, which is responsible for the quark parton model and (b), (c) QCD diagrams of first order 
in as that give contributions to F2 which depend on the quark and gluon content respectively of the 
proton. 

Effectively, the term resums the leading InQ 2 terms. That is the (ci^mQ 2 )™ contri- 
butions which, in an axial gauge, correspond to the sum of ladder diagrams (with n rungs) in 
which the transverse momenta of the emitted partons (gluons) are strongly ordered along the 
chain (i.e. Q 2 3> k\ T 3> ... ^> k\ T in the example shown in Fig. 2). The NLO contribution 
corresponds to the case when a pair of momenta are comparable ~ &j+it and we lose a 
power of \nQ 2 . That is sums up the ag\n n ~ 1 Q 2 contributions. When truncating the 
power series at a given power of as, say a™, the renormalisation of the fi introduces a scheme 
dependence of 0(a s l+1 ). Traditionally the MS scheme is used. 



The parton densities describe not only deep-inelastic scattering, but all hard scattering pro- 
cesses with incoming nucleons. As we have noted, the densities fi(x, Q 2 ) have to be determined 
by experiment at some scale Q 2 = Qq- The basic procedure is to parametrize the x dependence 
at some low Q%, but where perturbative QCD should be applicable, and then to evolve up in 
Q 2 using the NLO DGLAP equations to determine fi(x,Q 2 ) at all the values of x, Q 2 of the 
data. The input parameters are then determined by a global fit to the data. 

To be specific the 1994/5 MRS g |] and CTEQ J§ analyses took Q 2 = 4 GeV 2 for the 
input scale. We describe the MRS analyses. Similar results are obtained by CTEQ. The starting 
distributions are taken to be of the form 



for i = g, w val , d val and the (total) quark sea S. In practice not all of the parameters 
Ai, Xi, /3i, 7j, 5i are free. Three of the are determined by the flavour and momentum 
sum rules. Moreover we have some idea of the values of the j3i and A, from spectator counting 



2. Global analyses 



x fi{x,Ql) = Ai x~^{l-x)^ (l + ^y/x + Si x) 



(5) 
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Figure 2: DGLAP evolution, which at LO sums the (aslnQ 2 ) n contributions, corresponds for P gg to 
the sum of such ladder diagrams with the transverse momenta of the emitted gluons strongly ordered 
along the chain {Q 2 3> k 2 T 3> ... 3> k 2 T ). On the other hand for the BFKL In 1/x summation (which is 
discussed in section 4) the ladder must be regarded as an effective ladder diagram incorporating many 
different contributions. In this case the (a^ In l/x) n contribution comes from the strongly ordered 
configuration x -C £ n <C ... <C £1 but with the gluon hr values unordered. 



rules and Regge expectations respectively. The QCD coupling is also a free parameter. It is 
determined primarily by the scaling violations observed in the high precision BCDMS _F^ p,Md 
data in the region 0.35 < x < 0.55. 

The flavour structure of the quark sea S = 2(u + d + s + c + ...) is determined by data. The 
CCFR dimuon production data imply that the strange sea is suppressed by 0.5 relative to the 
u and the d sea distributions at Q 2 = 4 GeV 2 . The difference d — u is arranged to be compatible 
with the observed NA51 || asymmetry in Drell-Yan production in pp and pn collisions. The 
input charm sea is determined by EMC deep-inelastic data for F 2 C 0. We assume c = for 
Q 2 < m 2 and for higher Q 2 we generate c(x, Q 2 ) by massless evolution. The data imply m 2 = 
2.7 GeV 2 . After evolution to Q 2 = 4 GeV 2 we find that the charm sea, to a good approximation, 
satisfies 2c = 0.02S. The description of the EMC charm data by MRS(A) partons is shown in 
Fig. 3. Clearly this is an approximate way to treat m c ^ effects. At this meeting De Roeck 
H presented the first preliminary measurements of F£ by the HI collaboration. To get some 
idea of the future impact of these data, estimates of the preliminary measurements of F% at Q 2 
= 13, 23 and 50 GeV 2 have been superimposed on the x = 0.002 curve in Fig. 3. Although 
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the present parton treatment of charm appears to be satisfactory, it is clear that future more 
precise data will be invaluable in the investigations of the proper treatment of m c ^ effects. 
In summary, the data imply that at the input scale, Qq = 4 GeV 2 , the charm sea carries about 
0.4% of the proton's momentum, as compared to nearly 4% by the strange sea, 6% by the up 
sea and 9% by the down sea. 




5 10 50 100 
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Figure 3: The description of the EMC measurements of F% with x > 0.024 by the MRS (A) partons 
The preliminary measurements of F% by the HI collaboration, presented for the first time at this 
Conference J|], have been used to superimpose three data points to be compared with the x = 0.002 
curve. 

The wide range of data used in the global fits is shown in the top part of Table 1, together 
with an indication of the most important constraints that they impose on particular partons. 
Fig. 4 shows the parton distributions at Q 2 = 20 GeV 2 corresponding to the 1994 and 1995 
sets of MRS partons 0, §], which were obtained from global fits to these data. The differences 
u — u = w va i and d — d = d va \ show the valence quark structures around x ~ 0.1. The dominance 
of the gluon for x < 0.01 is also evident. Note that the gluon is suppressed on the figure by 
a factor of 10. The general conclusion is that the partons are well determined for x > 0.02, 
where data for a wide range of processes exist, except possibly the gluon which in this region 
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Table 1: The experimental data (in the top part of table) used to determine parton distributions in 
the global analyses. The last column gives an indication of the main type of constraint imposed by a 
particular set of data. The processes in the bottom part of the table are discussed in Section 3. 



Process and Experiment 



Leading-order 

subprocess Parton and 0:5 determination 



DIS (/jlN -> fiX) 7*g -> q 

BCDMS, NMC, E665 
Ff , F% n 

DIS (uN -> fiX) W*q -> q' 
CCFR (CDHSW) 

/xN — > ccX 7*0 — > c 

F 2 C , EMC 

i/iV -> X W*s -> c 

CCFR /i+ 

DIS (HERA) 7 *g -> g 

F 2 ep (H1,ZEUS) 

pp -> -yX WA70 (UA6, g# -> 75 

E706, R806, UA2, CDF) 

pN — ► n + n~X qq — > 7* 

E605 

pp,pn — > p + p X tttt, c?J— > 7* 

NA51 ud,d«->7* 

pp — >■ W ± asym u J — > IU + 

CDF cfe -> W" 



Four structure functions - 

u + u, d + d 

u + d, s (assumed = s) 
but only / xg(x)dx ~ 0.5 
[u — d is not determined] 
as (x ~ 0.4 data) 

c « 0.1s at Q 2 = 4 GeV 2 
s ~ \u (or |d) 
A 9 , Ag, as 

(xg ~ x~ Xg ,xq ~ x~ As ) 
(7(2; 0.4) 

q = ...(l-xf 3 

(u- d) at x = 0.18 

slope of u/d at rr 0.05 



pp -> jets #g, qq #(0.005 < x < 0.1) 

CDF, DO ^2 jets q(x ~ 0.2) 

as (Ex ~ 100 GeV) 

7*p -> dijets 7*3 -> gg #(0.005 < x < 0.1), a 5 

HI, ZEUS 7*g -> #g 

7*p -> 7*3 -> (cc)^ #(?) 
EMC, HERA 

7P — > J/i/>P cc — > cc g(x ~ 10~ 3 ) 

HI, ZEUS via gg exch. 
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Figure 4: The 1994 @ and 1995 § sets of MRS partons at Q 2 = 20 GeV 2 . For clarity the gluon 
distribution (divided by a factor of 10) is only shown for x < 0.005. 

is mainly constrained by prompt photon data. The gluon is clearly the crucial parton in the 
small x domain. It is the subject of the next two sections. 



3. Determination of the gluon 

The gluon only contributes at leading order in prompt photon production among all the 
processes fitted in the global analyses, see Table 1. Its distribution is therefore not so well 
determined as those of the quarks. The constraints on the gluon come mainly from (i) the 
momentum sum rule, (ii) prompt photon production, and (iii) the scaling violations of F 2 . 
Scaling violations impose the tightest constraint in the small x region, where the gluon is the 
dominant parton. Then 

dF 
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Figure 5: The x intervals in which the gluon may be constrained by various sets of data. Also shown 
are the gluons at 20 GeV 2 from ref. Q. 



where the convolution leads to the gluon being sampled at a higher value of x than that at 
which the violation is measured. Roughly speaking an observed violation at x measures as{Q 2 ) 
g(2x). In this way the HERA measurements of the scaling violations of F2 at small x have 
considerably improved our knowledge of the gluon. 

These and other potential determinations of the gluon are summarised in Fig. 5, with an 
indication of the relevant x ranges. We discuss the determinations in turn below. 



3.1 F2 scaling violations at HERA 

In the early NLO global fits not all the parameters in 

xg = A g x~ x » (1 - xf 9 (1 + 7 5 y/x + 5 g x) (7) 

were used. In particular, the data did not determine the small x behaviour of the gluon. For 
example, 7 g was set to zero, and since in the perturbative region the gluon drives the sea, via 
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g — > gg, it was assumed that A 9 = As at the input scale Qq. With the advent of the HERA 
measurements of F%, and their improvement year-by-year, the gluon has become better and 
better determined in the small x region. The improvement is reflected in Table 2 in the step- 
by-step release of the parameters, X g and 7 9 , which most affect the small x behaviour of the 
gluon. The G set of partons allowed As 7^ X g for the first time, but it lead to only a marginal 
improvement with respect to the set A' in which X s was set equal \„ @. The new HERA data 



lT| |12|| exclude the G set of partons and it is interesting to note that the new parton set Rl 



TOf is similar to the A' set of partons. 



Table 2: The exponents Aj of xg ~ x 
sequence of MRS analyses 



Xg and xS 



Xs at the input scale Ql = 4 GeV 2 for a 



H|, [|, 10] which include more and more precise HERA data as they 
become available each year. The latest fit [1C] has input scale Qq = 1 GeV 2 but the exponents are 
also given after evolution up to Q 2 = 4 GeV 2 . The values of the parameters Aj are strongly correlated 
with the values of 7$ (see the discussion in ref . M ) . 



MRS fit X g A 



1993 


D 








fixed 


(l 9 = o) 




D_ 


0.5 


0.5 


fixed 


(7, = 0) 


f994 


A 


0.3 


0.3 


free 


{l 9 = o) 


f995 


A' 


0.17 = 


0.17 


free 






G 


0.31 


0.07 


free 




f996 


Rl 


(0.17) 


(0.18) 


at Q 2 


= 4 GeV 2 






-0.55 


0.12 


at Qq 


= 1 GeV 2 



Traditionally the MRS and CTEQ analyses have fitted to data with Q 2 > 5 GeV 2 . In the 
GRV approach JI3| , [TJJ valence-like forms of the parton distributions are taken at a low input 



scale Qq = 0.34 GeV 2 . The original hope of this "dynamical" model was that input valence 
quarks would suffice and that the gluon and sea distributions would be generated radiatively. 
However, a sizeable valence gluon and a valence sea distribution are also required at the input 
scale in order to describe prompt photon and NMC deep-inelastic data respectively, which as 
a consequence introduces more phenomenological parameters into the GRV model. The GRV 
partons were found to give a good description of the HERA data down to unexpectedly low 
values of Q 2 , namely Q 2 ~ 1.5 GeV 2 , although with the precision of the latest data there is 
some discrepancy at the lowest values of x, see Fig. 6. Motivated by the general success of the 
GRV (DGLAP-based) predictions, the latest MRS analysis |TU| uses a lower input scale, Qq = 1 
GeV 2 , and fits to data Q 2 > 1.5 GeV 2 — the resulting description of the new HERA data at 
the lowest values of x is shown in Fig. 6. The continuous and dashed curves correspond to 
setting as(Mf ) = 0.113 and 0.120 respectively. The first choice of as is the value determined 



S 



by the scaling violations of the fixed-target deep-inelastic data, in particular the BCDMS F^'^ 
measurements in the interval 0.35 < x < 0.55 |]15|. The second choice is preferred by LEP data 
||16|| , and also marginally by the HERA measurements of _F 2 ep . The two new MRS fits have both 
X g and Xs as free parameters. To see the extent to which they can be determined independently, 
fits are also performed with X g = Xs- The gluon distributions of these 4 fits at Q 2 = 5 GeV 2 
are shown in Fig. 7. Also shown in this plot is a representative spread of the gluons that were 
available in 1995. The large difference between the A' (with As = X g ) and G (with Xs ^ X g ) 
gluons is not present in the new fits — they all cluster about A'. The new HERA data appear 
to have significantly pinned down the gluon. 



New NLO fits and. GRV compared, with new HERA data 




Figure 6: Recent HERA measurements [|ll], ^] of the proton structure function F2 versus InQ 2 at 
the lowest values of x. The continuous and dashed curves are the description obtained in a new 
preliminary MRS global analysis with the QCD coupling taken to be such that as(M^) = 0.113 and 
0.120 respectively [10|. The dotted curve is the prediction obtained from the GRV partons [13]. 
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Gluon densities 




Figure 7: The gluon distribution at Q 2 = 5 GeV 2 . The four continuous curves show a representative 
set of the gluons that were available in 1995. The four broken curves correspond to gluons obtained 
[ pj] ] in global analyses which include the latest HERA data — the new solutions cluster around the 
MRS(A') gluon. 

3.2 Prompt photon production 

The processes pp — > 7X and pp — > 7X have long been regarded as a classic way to determine 
the gluon. The perturbative QCD formulae are now known to NLO, including the fragmentation 
(or Bremsstrahlung) contribution [fT? ]. Moreover, the experiments (WA70, UA6, E706, R806, 



UA2, CDF) cover the entire x interval from 0.6 down to 0.01. The situation therefore appears 
promising. However, there is a pattern of deviation in the shape of the pt dependence. The data 
are steeper in p? than the QCD predictions. Neither changes of scale nor the introduction of 
fragmentation effects can resolve the discrepancy! since the various experiments probe different 



ranges of x ~ xt = %Pt/ 'y/s. On the other hand it has been shown [ll| that the discrepancy 



2 Vogelsang and Vogt |Q have demonstrated that these effects can improve the description of a single 
experiment. However, experiments at different y/s reproduce a similar pattern, but in different x intervals. 
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can be removed by a broadening of the transverse momenta of the initial state partons (due to 
multigluon emission) which increases with the energy y/s. A similar effect has been quantified 
in Drell-Yan production, but here, so far, the broadening is accounted for phenomenologically. 
Until the multigluon effects are calculated in QCD it is not possible to use the prompt photon 
data (especially those at higher energy) to pin down the gluon. The most reliable determination 
comes from the lower energy pp — > 7X data of the WA70 collaboration, where the broadening 
is much less!, but even here there is an ambiguity of some ± 25% in the value of the gluon. 



3.3 Jet production at Fermilab 

Dijet production in pp collisions can also, in principle, probe the small x behaviour of the 
gluon []20], ^T[. For example, if the two jets are produced with equal transverse momentum pt 
but both very forward with pseudorapidity r\ ^> 1 then x\ ~ 1 and x<i ~ {2pt/ \fs) exp(— 77) <C 1. 
Detailed NLO calculations [^TJ show that at y'i = 1.8 TeV the gluon can be probed in this 
way in the range 0.005 < x g < 0.05. However, at present the systematic errors are too large to 
allow any definite conclusion to be drawn. 

The single jet inclusive cross section for jet transverse energies Er ~ 100 GeV is dependent 
on the gluon via the gg,gq,gq initiated subprocesses. The gluon g(x,Q 2 ) is sampled at x ~ 
2E T /y/s ~ 0.1 (and Q 2 ~ E?f) for centrally produced jets. However, the E T spectrum gives 
more information on the running of as{E^) than on the gluon. In fact if we take the information 
on the gluon from the scaling violations of F2 at HERA, then the jet Et spectrum gives a 
sensitive measure of as |22]. The steeper the spectrum the larger the prediction for as- There 



are indications from the medium E? CDF (and also from the preliminary DO) jet data that 
the observed spectra favour as(Mf ) in the region 0.116 to 0.120 [|2|, |T0|. Again the systematic 
error is the limiting factor. 



3.4 Dijet production at HERA 

The observation of dijets in deep inelastic scattering at HERA offers, in some respects, 
similar possibilities to jet production at Fermilab. Again within a single experiment it is possible 
to observe the running of as(fcy). At HERA the LO subprocesses are the QCD Compton process 
11 ~^ gii an d, relevant for the gluon, the •yg — > qq fusion reaction. The NLO contributions 
are known and the scheme dependence has just been quantified. Indeed Mirkes and Zeppenfeld 



23[ have presented to this Conference a detailed study of the jet algorithms and conclude that 
the cone or kx schemes are favoured and lead to less scale dependence than the other schemes. 
The clean identification and kinematic measurement of the jets is the experimental challenge. 



3 Also pp — > 7X has the advantage that the dominant LO subprocess is gq — > jq, unlike pp — > yX where 
qq — * jg gives a comparable contribution. 
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3.5 Inelastic Jjip photoproduction 

It has long been advocated that inelastic J/ip photoproduction at HERA may serve as a 
measure of the gluon — see, for example, ref. [24j which considers the colour-singlet model [25 
for the process at LO accuracy. Recently the NLO contributions have been calculated [Ef 



A detailed study of the spectra in the high energy range at HERA shows that the perturbative 
calculation is not well-behaved in the limit px — > 0, where pt is the transverse momentum of the 
J/ip. No reliable prediction can be made in this singular boundary region without resummation 
of large logarithmic corrections caused by multigluon emission. If the small pt region is excluded 
from the analysis, the NLO result accounts for the energy dependence of the cross section and 
for the overall normalization, see Fig. 8 [27|. However, since the average momentum fraction 



of the partons is shifted to larger values when excluding the small-pr region, the sensitivity of 
the prediction to the small-x behaviour of the gluon distribution is not very distinctive. 
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Figure 8: Total cross section for inelastic J ftp photoproduction as a function of the photon-proton 
energy for different parametrizations of the parton distribution in the proton. Experimental data from 
pqj. The figure is from [27]. 



3.6 Diffractive J/tp production at HERA 



Diffractive J/ip photoproduction appears to offer a more promising way to distinguish be- 
tween the gluon distributions. Since this is essentially an elastic process the cross section is 
a measure of the square of the gluon density. To leading order the cross section is given by 



29, 30 



J/ipp) 



48a 
12 



— =8 — [xg{x, Q )] 



(8) 



with Q = jM 2 and x = MZ/W , where W is the 7p cm. energy. In a recent study pi 
corrections to this formula have been calculated and comparisons with HERA data made, see 
Fig. 9. It was emphasized that the W dependence, rather than the normalisation, was the more 
reliable discriminator between the gluons. The power of the method is evident from Fig. 9, 
which appears to favour the MRS (A') gluon. Further phenomenological studies of this process 
can be found in refs. 1132 



a 



> 

T 



★ HI 94 Prelim. 




O ZEUS 94 Prelim. 




■ ZEUS 93 




• E401 




r A E516 


,4 W 


/V 


GRV(94) 




MRS(A') 


: f 


MRS(G) 



W (GeV) 



io 2 



Figure 9: The measurements of the cross section for diffractive J/ip photoproduction compared with 
the full perturbative QCD prediction obtained from three 1994/5 sets of partons. The figure is taken 
from 1 31]. 



4. The gluon at small x 

We have seen that the gluon is by far the dominant parton in the small x regime. Indeed 
in the perturbative region it drives the entire partonic structure of the proton via the g — >■ gg 
and g — >• qq transitions. 

So far our description of the data, including the small x HERA data down to Q 2 = 1.5 
GeV 2 , has been based on the DGLAP resummation of LO and NLO InQ 2 terms. In fact the 
rise of F2 with decreasing x that is observed at HERA appears to be well described by the 
simplest approximation for the small x behaviour of the gluon. If in the DGLAP evolution for 
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gluon, we take the splitting function equal to its small x limit, P gg ~ (3as/ir)/x, then it follows 
that 



xg(x, Q 2 ) ~ xg(x, Ql) exp 2 



55 In (*■) In ( i 

25 U / Vx 



(9) 



where t = ln(Q 2 /A 2 ), modulo slowly varying logarithms, provided the input xg(x, Qq) is not 
singular. That is, in this double logarithm approximation, xg increases faster than any power 
of ln(l/x), but slower than a power of (1/x). This behaviour feeds through into F 2 and gives 
an excellent description of the data, as emphasized by Ball and Forte [Q; the steepness of the 
rise in F 2 can be tuned to the data by adjusting the evolution length Q 2 /Qf ) . 

Does the success of the DGLAP description of the HERA data indicate the dominance 
of the InQ 2 resummations and the absence of higher twists? Such a conclusion would be 
premature. At small x, x < 1CT 3 , we have, so far, only one type of data (F 2 ep ) and there is 
freedom in the description, particularly as we have to supply the non-perturbative input at 
some scale Qq. Clearly at sufficiently small x the (NLO) DGLAP evolution will break down. 
When atslnl/x ~ 1 we have to also resum a^lnl/x contributions (or, to be more precise, a 
whole series of ln m Q 2 ln n \jx terms). At LO, the (aglnl/x)™ resummation is accomplished 
by the BFKL equation ^5] . In a physical gauge, the (a^lnl/x)™ term corresponds to an n- 



rung effective ladder diagram (Fig. 2) in which the soft gluon emissions are strongly-ordered 
in longitudinal momenta, but in which the transverse momenta are no longer ordered. Due 
to the latter fact we have to introduce the gluon distribution f(x,k^) unintegrated over k?, 
and anticipate a diffusion in Ink 2 -, as we proceed along the gluon chain. The relation to the 
conventional gluon is given by 

rQ 2 rlk 2 

xg(x,Q 2 ) = / -ff(x,k 2 T ). (10) 
In principle, in the small x domain the unintegrated distributions are the universal parton 



distributions which link process to process, via the kr factorization theorem [|36|. An example 
of the theorem is given below in (|T3| ) . 

For fixed as the x — > behaviour of the BFKL solution can be written in analytic form. 
Keeping only essential factors it behaves as 

where A^, the famous BFKL intercept, is given by Xl = (3as/7r)4 In 2. Thus we have a x~ Xl 
power-like growth accompanied by a diffusion in \nk^. If a physically reasonable prescription 
for the running of as is assumed then the BFKL equation may be solved numerically to yield 



37] a form 

f{x,k 2 T ) ~ C(k 2 T )x- x (12) 

where A ~ 0.5 is less sensitive to the phenomenological treatment of the infrared region, than 
the normalisation C. The BFKL prediction for the structure function F2 is obtained using the 
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krp factorization theorem 36 



f 2 (*,q 2 ) = j ^/(x',4)^g,4,Q 2 ), (13) 

where the integration is over the kinematic variables (x', k^) of the virtual gluon coupling to 
the quark box in Fig. 2. F 2 9 is the off-shell gluon structure function which at LO is given by 
the quark box (and crossed-box) contributions to photon-gluon fusion. The BFKL approach is 
also found to give a satisfactory description of F2 for small x. However, there are, at present, 
limitations to the "prediction". We comment on the ambiguities in the BFKL calculation of 
F 2 below. 

(i) Due to the diffusion of f(x, fcfi) in In k\ there is a significant contribution from the infrared 
k\ region which is beyond the scope of perturbative QCD and which has to be included 
using physically motivated phenomenological forms. This leads to an uncertainty in the 
overall normalization of F 2 , but much less in the x dependence. A physically reasonable 
treatment of the infrared region is found to give the experimental normalization. In a 
sense this is equivalent to providing the non-perturbative input for DGLAP, but here the 
x~ x behaviour at small x is prescribed. 

(ii) The BFKL equation only resums the LO In 1/x terms. The NLO contributions are needed 
for a stable prediction. Sub-leading effects, which to a large extent embrace both energy- 



momentum conservation and angular ordering, have been shown [^, [38| to significantly 
reduce the value of the exponent A. 

(iii) An underlying soft Pomeron contribution has to be included in the small x region, deter- 
mined by the extrapolation of the observed values of F 2 at large x. Again this has the 
effect of reducing the value of A apparent in F 2 ~ x~ x . 

(iv) Shadowing corrections to the BFKL equation will eventually, as x decreases, suppress 
the x~ x growth. Although they have not yet been fully formulated, the evidence from 
the observed ratio of diffractive to non- diffract ive deep-inelastic events, and from the 
persistent rise of F 2 at very low Q 2 ~ 2 GeV 2 , indicates that shadowing effects are at 
most 10% in the HERA regime. 

(v) We need further studies of a unified approach which incorporates, on a sound theoretical 
footing, both the BFKL and DGLAP resummations. 

From the above discussion it is clear that there are many issues to be resolved. We see that 
it will not be easy to quantify the importance of the In 1/x BFKL- type contributions by using 
the effective A dependence of the measured values of F 2 at small x, that is F 2 ~ x _Acfi . 

There has recently been much activity |39|] based on expanding the anomalous dimensions 



jij^agjU)) in terms of as and the moment variable u. For instance, for the gluon anomalous 
dimension, (LO) DGLAP resummation amounts to summing the terms 

a s 

Igg = d x h d 2 a s + d 3 a s uJ + . . . , (14) 

UJ 
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whereas BFKL resums a different subset of terms 

7 33 = &i — + &4 4 + • • " ( 15 ) 

UJ UJ^ 

with the uj~ n term corresponding to a x _1 log n_1 l/x contribution to P gg (except for n = 1 
where P gg ~ 1/x, see (|T6|)). Both the above expansions start with the same "double logarithm" 
term 

rl rl ' diets \ as 



7 S9 = J dxx u P gg ~ J dx x^ y j = d\ — , (16) 

with b± = d\ which leads to the behaviour displayed in (M). Some of the first few coefficients of 



the BFKL expansion (|15|) vanish, namely 62 = &3 = ^5 = [|0|. For this reason much of the 



rise of F 2 with decreasing x is attributed to the LO expansion of 7 gg 

Iqg = a S |Ci + C 2 — + C 3 — + ...J, (17) 

which contributes to F2 via <9i<2 /<9 In Q 2 = P qg ® g. All the coefficients q are non- vanishing 



(unlike those for 7 gg ), positive definite and large [41]. 

In principle, this approach appears to offer the attractive possibility of quantifying the 
importance of In 1/x effects by studying DGLAP-type evolution with anomalous dimensions 
(and coefficient functions) which incorporate the (as/uj) n terms. However, it has been pointed 
out |^2| that such a procedure masks the true dependence on contributions from the infrared 
region. Due to the diffusion in lnk^ at small x, an (as/uj) n contribution, which in DGLAP is 
assigned to the local point Q 2 , actually samples the region of (the logarithm of) virtuality 

InQ 2 ± (A(x))" (18) 

where A increases approximately as (In 1/x) 2. Due to the large numerical coefficient, B = 
(3as/7r) 56^(3) (with £(3) = 1.202), in the diffusion term in (TlT|), there is considerable implicit 
penetration into the infrared region. For example it is found that the (as/uj) 4 term samples 
virtualities down to Q 2 /100 ||. Thus [t 

seems that there is no alternative but to work with 
the unintegrated gluon distribution and the kr factorization theorem and to study the effects 
of contributions from the infrared kx region explicitly. 

The observable F2 is too inclusive to show all the characteristics of the small x properties 
of the unintegrated gluon f(x, k^)- In particular the lnA;^ diffusion pattern is integrated over. 
For this reason other observables which measure properties of the final state in deep inelastic 
scattering have been advocated as better indicators of In 1/x resummation effects (see, for 
example, the reviews listed in ref. 



5. Conclusions 

The universal parton distributions of the proton are well determined by a wide range of 
data in the region x > 0.02. The exception is the gluon, which although constrained, still has 
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some residual ambiguity. However, the scaling violations observed in the new, more precise 
HERA measurements of F 2 have pinned down the gluon in the small x region (x ~ 10~ 3 ). 
One consequence is that the GRV model which gave such an excellent description down to 
Q 2 ~ 1 GeV 2 , shows some systematic discrepancy. The prediction is above the new HERA 
measurements at the lower values of x, see Fig. 6. Indeed the previous spread of possible gluon 
behaviour at small x (represented by the MRS (A', G), CTEQ3 and GRV curves in Fig. 7) has 
been narrowed in global analyses 0] incorporating the new data to give gluons similar to that 
of the MRS(A') set. Motivated by the success of the GRV approach, the latest global analysis 
[TOf uses a lower input scale, Qq = 1 GeV 2 . 

The measurements of diffractive J/ip photoproduction at HERA were seen to also favour 
the gluon of the MRS (A') set of partons, see Fig. 9. We briefly reviewed other ways in which the 
gluon may be measured. Jet production at Fermilab and at HERA offer not only a constraint 
on the gluon, but also provide a sensitive measure of the running of a$- 

Finally, we briefly discussed the perturbative QCD expectations for the behaviour of the 
gluon in the small x region. We emphasized the importance of resumming the (a^lnl/x) 71 
terms when x is sufficiently small so that a^lnl/x ~ 1. We highlighted the problems of 
incorporating these effects in the description of F 2 and the necessity to use the unintegrated 
gluon distribution together with the factorization theorem. The dramatic improvement in 
the experimental measurements in the small x domain serves as a challenge to provide a deeper 
theoretical understanding of this fascinating frontier of QCD. 
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